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Abstract 

A very promising spin physics programme will be soon on the way at the BNL Relativis¬ 
tic Heavy Ion Collider (RHIC). By studying the spin asymmetries for various processes 
(single photon, single-jet and W ± production), we will compare the different predictions 
obtained using some sets of polarized parton distributions, available in the recent liter¬ 
ature. We will put some emphasize on the analysis of the anticipated errors, given the 
event rates expected from this high luminosity new machine and the current acceptance 
for the detector systems at RHIC. 
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1 Introduction 


Serious progress have been made over the last ten years or so, in our understanding of 
the spin structure of the nucleon. This is partly due to the discovery in 1987, by the EMC 
experiment [lj], of the violation of the Ellis-Jaffe sum rule |2j. This result was obtained by 
measuring the proton polarized structure function g{(x, Q 2 ), in polarized deep-inelastic 
scattering with a polarized amonia target ( NH 3 ). It led to an unexpected situation, that 
is, only a small fraction of the proton spin is carried by quarks, which is at variance with 
the naive quark-parton model. The EMC result has motivated several experiments, in 
three different laboratories CERN , DESY and SLAC J3J, to undertake some systematic 
determinations of gi’ n (x, Q 2 ), for proton and neutron by means of measurements using 
a polarized lepton beam and several different polarized targets (hydrogen, deuterium, 
helium-3). The main goal of these experimental programmes was to confirm the EMC 
result on gi(x,Q 2 ) and to measure the neutron structure function g™(x, Q 2 ), in order to 
test the fundamental Bjorken sum rule [£|. It turns out that these fixed target experiments 
allow to cover only a rather limited kinematic domain that is, 0.005 < x < 0.7, with the 
corresponding average Q 2 , < Q 2 > between 2 GeV 2 and 10 GeV 2 . Needless to recall 
that the low x region, which is not easily accessible and therefore loosely determined, is 
very crucial for testing accurately these sum rules. In addition it is not possible to make 
a flavor separation to isolate the contribution of each quark to the nucleon spin. 

On the theoretical side, the surprising EMC result has also stimulated an intense 
activity, in particular its possible interpretation in terms of the so-called axial anomaly 
which involves the gluon polarization A G(x, Q 2 ) BS0. This is an interesting physical 
quantity in itself, but it can only be determined indirectly in polarized deep-inelastic 
scattering. More precisely the contribution of the gluon polarization to the Ellis-Jaffe 
sum rules is ambiguous and, as we will see below, much cleaner information will be 
extracted at the RHIC facility from collisions with polarized proton beams. Let us also 
mention that several important results in higher order perturbative QCD corrections have 
been obtained recently a as- which will allow, in the future, very stringent tests of the 
QCD spin sector, so far essentially unexplored. 

In order to describe the most recent experimental results on g\’ n,d (x , Q 2 ), several phe¬ 
nomenological models have been proposed and correspondingly, one can find in the recent 
literature various sets of parametrization for the polarized parton distributions. In spite 
of the constant progress in the accuracy of the data, the extraction of these distributions, 
especially for sea quarks and gluons, remains fairly ambiguous and controversial. We 
will select some of them and recall their main features in Section 2, where we also give 
some basic kinematics and the definitions of the spin observables we will study below. In 
Section 3 we will consider direct photon production. In terms of perturbative QCD, it is 
a simple process whose double helicity asymmetry is particularly sensitive to A G(x,Q 2 ), 
as we will illustrate by means of a useful approximate expression. Section 4 is devoted to 
single-jet production which has a higher complexity than the previous process but whose 
generous cross section will allow a very accurate determination of the double helicity 
asymmetry. We will consider in Section 5 the production of W ± gauge bosons and parity 
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violating asymmetries with either one or two polarized beams. Their measurements will 
provide a good calibration of the quark and antiquark polarizations and also the u and 
d flavor separation. We will comment on the possible benefit one can get by measuring 
them simultaneously in proton-proton and proton-neutron collisions. We will give our 
concluding remarks in Section 6. 


2 Spin observables and various parametrizations 


Fundamental interactions at short distances which are probed in high energy hadronic 
collisions, involve hard scattering of quarks, antiquarks and gluons. Let us consider the 
general hadronic reaction 

cl T b —> c f A (1) 

where c, in the cases we will consider below, is either a photon, a W ± or a single-jet. In 
the hard scattering kinematic region, the cross section describing (1) reads in the QCD 
parton model, provided factorization holds, as 


da(a + b —> c + X ) = ^ 


1 + 8ij 


dx a dxb 


fl a \xa,Q' 2 )fi°\ x b,Q' 2 )dd v + (i^j) 


p(6) 


( 2 ) 


The summation runs over all contributing parton configurations, the /(x, Q 2 )’s are the 
parton distributions, directly extracted from deep-inelastic scattering for quarks and anti¬ 
quarks and indirectly for gluons, day,- is the cross section for the interaction of two partons 
i and j which can be calculated perturbatively, some of which , at the lowest order, are 
given in ref.||. If we consider the reaction (1) with both initial hadrons, a and b longi¬ 
tudinally polarized, one useful observable is the double helicity asymmetry All defined 
as 

da a{+)b{+) - da a(+)fe( _) 

All ~ l -rj-’ 

«° r a(+)6(+) + «0 r a(+)6(-) 


when we assume parity conservation, i.e. dcr a (A)fe(A') = d<7 a (_A)&(_A')- Its explicit expression, 
assuming factorization, is given by 


A 


LL 


da = J2 

ij 


1 + Sij 


dx n dxh 


A// a) (x a , Q 2 )Aff ) (x b , Q 2 )al L dd 13 + (*<-► j) 


( 4 ) 


where da is given by eq.(2) and al L denotes the corresponding subprocess double asym¬ 
metry for initial partons i and j. The A/’s are defined as 

A/(x,Q 2 ) = f+(x,Q 2 ) + f-(x,Q 2 ) (5) 

where f± are the parton distributions in a polarized hadron with helicity either parallel (+) 
or antiparallel (-) to the parent hadron helicity. Recall that the unpolarized distributions 
are / = /+ + /_ and A/ measures how much the parton / ’’remembers” the parent 
hadron helicity. If the subprocess involves parity violating interactions, one can consider 
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another interesting observable which requires only one initial hadron polarized, that is 
the single hclicity asymmetry Al , defined as 


Ai, = 


da a{+) 
d(j a (—) A d<r a ( + ) 


( 6 ) 


In addition, if both a and b are polarized one can also have two double hclicity parity 
violating asymmetries defined as 


A PV - 

n LL — 


°(+)K+) 


d(?g(-)b {-) ~ da, _ 

dcr a (-)b(-) A dCa(+)6(+) 


and 


All = 


da n( — 


«(-)&(+) 


da, 


'«(+)&(-} 


do"a(—)6(+) A dcT a (-|_)b(—) 


( 7 ) 


which can be simply related to Al as we will see in section 5, where we will give their 
explicit expressions for W ± production. This completes the presentation of the spin 
observables and let us now turn to a short discussion on the selected sets of polarized 
parton distributions. 

We first consider the set of polarized parton densities which has been proposed in 
ref.@. The quark and antiquark densities are parametrized in terms of Fermi-Dirac dis¬ 
tributions with only nine free parameters. Simple assumptions allow to relate unpolarized 
and polarized quark distributions and the free parameters are determined from the very 
accurate data on deep inelastic neutrino scattering at Q 2 = 3 GeV 2 . For the gluon 
distributions one uses a Bose-Einstein expression which is given with no additional free 
parameters. This allows a straightforward DGLAP Q 2 evolution which leads to excel¬ 
lent predictions for F- 2 (x, Q 2 ) up to the x, Q 2 kinematic range accessible at HERA. These 
quark densities give also a fair description of the low Q 2 structure functions g± ,n (x,Q 2 ). 
However, in this approach A G(x, Q 2 ) at Q 2 = 3 GeV 2 is not known and two possibilities 
have been proposed, namely the soft and the hard gluon polarization which are shown in 
Fig.la. 

There is another standard method for obtaining the polarized parton distributions, 
which is more generally used in the literature. It is based on the direct analysis of 
gi’ n (x, Q 2 ), independently of the unpolarized structure functions. A very simple 
parametrization for the quark densities, of the type x a (l — x )A has been proposed in 
ref. [|10[. It was found that for a correct description of the data, in particular the very 


recent accurate neutron data [|TT[] from SLAC E154 at Q 2 = 5 GeV 2 , one needs a rather 
substantial A G(x,Q 2 ), having also a simple expression as above, which is displayed in 
Fig.la. 

In ref.|l^] the parametrization is based on more complicated expressions of the type 
x“(l — x)^(l + 7 \fx A px) in a similar way as what was used by Martin-Roberts-Stirling 
to describe the unpolarized quark distributions JR|. Due to our poor knowledge of 
A G(x,Q 2 ), they parametrize it in the same way and they explore various possibilities 
with three different choices of jc and Pg- They conclude that three scenarios A, B and C, 
shown in Fig.lb are equally acceptable and therefore AG(x, Q 2 ) is largely undetermined. 
However their analysis, which is prior to the SLAC E154 neutron data, does not include 


it. This is also the case for the parametrization proposed in ref.[14], where the polarized 
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distributions are constructed from the unpolarized ones, the so called GRV densities [|15 
by means of a simple multiplicative factor x a (l — x) b . Needless to recall that the free 
parameters are such as to respect the fundamental positivity constraints down to the low 
resolution scale of this approach namely Q 2 = 0.23 GeV 2 . In this case, one gets a gluon 
polarization, shown in Fig. lb, which is half way between the soft and the hard AG of 
ref- 0 • 

In addition to these four sets of polarized parton distributions we just briefly discussed, 
many other different choices have bee presented in the literature [|l^, [L7| . It appears that 


from the present available data, the valence quark distributions are fairly well determined 
and rather consistent with each other, unlike for the sea quarks (or antiquarks) and for 
gluons as we have seen above. 

Measurements of polarized deep inelastic lepton nucleon scattering yield direct infor¬ 
mation on the spin asymmetry 




N 


g?(x,Q 2 ) 

Ff(x,g 2 ) 


(■N = p,n,d ) , 


( 8 ) 


where F^(x,Q 2 ) is simply related to the unpolarized structure function F^(x,Q 2 ). All 
data for Q 2 > 1 GeV 2 are consistent with a flat Q 2 behaviour of A^(x, Q 2 ) which reflects 
the fact that the scaling violations of gi(x,Q 2 ) and F^{x,Q 2 ) are not so different and 
not distinguishable within experimental errors. The next-to-leading or leading order QCD 
evolutions leads to A±(x, Q 2 ) almost independent of Q 2 in accordance with experimental 
observation. To illustrate this fact we show in Fig.2 the results for the proton case 
Ai(x, Q 2 ) from ref.[[T3J. This interesting property will be used in the next section on 
direct photon production. 


3 Direct photon production 

The cross section for direct photon production on pp collisions at high px is considered 
as one of the cleanest probe of the unpolarized gluon distribution G(x, Q 2 ). This is partly 
due to the fact that the photon originates in the hard scattering subprocess and is detected 
without undergoing fragmentation. Moreover in pp collisions the quark-gluon Compton 
subprocess qG —» gy dominates largely and the quark-antiquark annihilation subprocess 
qq —> Gy can be neglected. Consequently the double helicity asymmetry A 1 LL (see eq.(4)), 
which involves in this case only one subprocess, becomes particularly simple to calculate 
and is expected to be strongly sensitive to the sign and magnitude of A G(x, Q 2 ). For the 
Compton subprocess, clll whose expression at the lowest order is given in ref.||, is always 
positive and such that auX6 cm = 90°) = 3/5 where 9 cm is the center of mass (c.m.) angle 
in the subprocess. 

Before we proceed, we would like to present an approximate expression of valid 
at pseudo-rapidity rj — 0, where the cross section reaches its largest values. To a good 
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approximation we can write 


A ll(v = 0 ) = 7 


3 AG(x Tl Q 2 ) 


5 G(x t ,Q 2 ) 


A i(xt, Q 2 ) , 


(9) 


where Xt = 2 Pt/\Z~s, Pt being the transverse momentum of the outgoing photon and 
y/s the c.m. energy of the reaction pp —> 'yX. In eq.(9) one should take Q 2 = p\- and 
A\ is the deep-inelastic spin asymmetry we discussed above (see eq.( 8 )). This expression 
shows clearly the dependence of AJ L on AG and, in order to test its validity, we have 
compared the exact calculation of A/ L , using eqs.(2) and (4) also with Q 2 = p^, to the 
above approximation. The results are depicted in Fig.3a for two sets of polarized parton 
distributions and show the usefulness of this simplified expression (9). A\ L increases with 
Pt due to the rapid growth of A T [ with x and follows the sign of AG and its magnitude 
which is larger in the case of BS-b than for GRSV (see Figs. la,b). 

In Fig. 3b we show a more complete comparison of the results we obtained for the 
Pt distribution of AC LL {j] = 0), using the different sets of polarized parton distributions 
we have discussed in the previous section, with a leading order Q 2 evolution. Actually 
we find that the smallest predictions correspond to the sets ref. 0 and ref. |L4| which 
have the smallest A G(x)/G(x). The predictions differ substantially at large Pt, which 
corresponds to the region, say around x = 0.4 or so, where the distributions A G(x) have 
rather different shapes. We have also indicated the expected statistical errors based on 
an integrated luminosity L = 800 pb ^ 1 at a/s = 500 GeV, for three months running time. 
We have evaluated the event rates in the pseudo-rapidity gap —1.5 < 7] < 1.5, assuming 
a detector efficiency of 100% and for a p/ acceptance A pj. = 5 GeV/c. We see that up to 
Pt = 50 GeV/c or so, A~l L w bl be determined with an error less than 5% which therefore 
will allow to distinguish between these different possible AG(x). For very large p/, the 
event rate drops too much to provide any sensitivity in the determination of AG(x). In 
Fig. 4 we display the pseudo-rapidity distribution of AA ll (i 7 ) at a fixed p/ value. It shows 
a systematic increase of AJ L for higher 77 , already observed in ref. R18 
due to the cross section fall off away from g = 0 . 


which is partly 


This smaller event rate leads also to 
larger statistical errors, but this measurement is worth performing because it is obviously 
a non-trivial QCD test. 

Let us now give a short discussion of these results and how they can be compared 
with other recent calculations. First, one may be concerned about higher order QCD 
corrections, to our dominant subprocess qg —> q 7 . The corresponding K-factors have 
been calculated [|19|, |20| 1 i and they always exceed 1, so they rescale up separately all 
spin-dependent cross sections. However this effect tends to be marginal on AJ jL which 
is a ratio of cross sections. Next we should consider the effect of a Q 2 evolution of the 
polarized parton distributions which were done here at the leading order (LO) while in 
ref. m they took the next-to-leading order (NLO) version. It is known that there is 
almost no difference for the valence (u and d) quark distributions, but this is not the case 


1 The results obtained in ref.| 
choice of a consistent scheme. 


which are slighty different from those of ref. p9|. are based on the 
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for the gluons. Actually at the NLO, the maximum of A G{x)/G{x) is shifted towards 
smaller x values and therefore the rise of AJ l is sharper. This is precisely what we see 
by comparing the GS sets in Fig.3b with Fig.4 of ref. [| 2 T[ , which also uses ref.||T 2 |]. Finally 


there is another aspect to be mentioned. At collider energies the outgoing photon is very 
often produced inside an hadronic jet, whose energy has to be restricted by isolation cuts 
in order to improve the detection accuracy. This is a standard procedure experimentally 
but not at all straightforward theoretically. It has been studied at the NLO with a Monte- 
Carlo method [|22| and it was found that the effects of isolation cuts are rather important 
for unpolarized cross sections, but not so much for AJ L and they decrease with increasing 
energy. Actually they become minute at t/s = 500 GeV and the predictions one finds in 
ref. [p 2 | are almost identical to those of ref. pl| where isolation cuts were ignored. 


4 Single jet production 

Inclusive jet production is also a physics area where one can learn a lot about par- 
ton densities and, considering the vast amount of unpolarized existing data, it has been 
regarded as an important QCD testing ground. Event rates are much larger than for 
prompt photon production, but there is a drawback because many subprocesses are in¬ 
volved, unlike in the previous case. In principle one should take into account gluon-gluon 
(GG), gluon-quark ( Gq ) and quark-quark ( qq ) scatterings. Although these subprocesses 
cross sections are not so much different, after convolution with the appropriate parton 
densities (see eq.( 2 )), they lead to very distinct contributions to the hadronic spin-average 
cross section. 

In the spin dependent case where no data exist so far, the same ingredients can be 
used to test earlier determinations of the polarized parton densities and to check the 
consistency of our picture of the nucleon spin. A first attempt in connection with RHIC 
was proposed in ref.JT 8 [, where various aspects of jet physics were discussed. Here we 
will restrict ourselves to the double helicity asymmetry Aj/' for single jet production 
and in order to clarify the interpretation of our results below, let us recall some simple 
dynamical features. In the very low p/ e< region, say p/ e< ~ 10 GeV/c or so, GG scattering 
dominates by far, but its contribution drops down very rapidly with increasing p/ et . In the 
medium range, say 20 GeV/c < p? 1 < 80 GeV/c or so, Gq scattering dominates and 
then decreases for large p? 1 , to be overcome by qq scattering. Of course these are rough 
qualitative considerations and accurate numerical estimates for the relative fractions of 
these different contributions depend strongly on the parton densities one uses. 

Let us first look at A^ l {i 7 = 0) and from the above discussion we see that one does not 
expect an approximate expression similar to eq.(9). However in the medium p range 
where Gq scattering dominates, A/^( 77 = 0) should have a trend similar to A/j/p = 0), 
with perhaps a magnitude reduced by a factor two, since about half of the jet cross section 
is due to GG and qq scatterings. This is what we see approximately in Fig. 5, where we 
present the numerical results for Aj^p = 0) at a/s = 500 GeV, which should be compared 
to Fig.3b. We have also indicated the statistical errors which are extremely small in this 
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case, because of the huge event rates. The p distribution of A 2 3 ^{p) for pT = 60 GeV/c is 
shown in Fig. 6 and it is either flat or decreasing for increasing rj. This is at variance with 
the prompt photon case (see Fig.4), because although the unpolarized jet cross section 
decreases away from p = 0 Q the numerator of A\j(r]) is sensitive to the behaviour of the 
different contributions GG, Gq and qq. We find that while Gq increases with p, both GG 
and qq are decreasing and therefore the precise behaviour of A^jijj) depends strongly on 
their relative sizes. The statistical errors, assuming a realistic A p = 0.2 per bin, will be 
small enough to allow the extraction of the correct trend from future data. Finally we 
recall that the predictions of ref. |2T| are somehow different from ours, which is certainly 
due to the use of different sets of polarized parton distributions. As already noted in 
the sets of refj 


ref. [21 


have different flavor decomposition and, as we will 


and ref. [12 

see in the next section, W ± production is certainly a very promising way to improve our 
knowledge on the u, d flavor separation. 


5 W ± production 


Let us first consider, for the reaction pp —> W ± X , the parity violating single helicity 
asymmetry Al defined in eq.(6). In the Standard Model, the W gauge boson is a purely 
left-handed object and this asymmetry reads simply for W ± production 


aT{v) = 


A u{x a , Mjy) d(x h , Myy) - {u^ d) 

u(x a , My?) d(xbi Myr) + (u d) 


( 10 ) 


assuming the proton a is polarized. Here we have x a = \Zre y , x& = \fre~ v and r = M^js. 
For W~ production the quark flavors are interchanged (u <-*■ d). The calculation of these 
asymmetries is therefore very simple and the results are presented in Figs.7a,b at y/s = 500 
GeV, for different sets of distributions. As first noticed in ref. |23|], the general trend of 
Al can be easily understood as follows : at y — 0 one has 


aT = 


1 Mu 


and 


u 


AT = 


( 11 ) 


evaluated at x — Mw/y/s = 0.164, for y = — 1 one has 


aW+ 

A l ~ 


Ad 

- j- and 

d 


evaluated at x — 0.059 and for y = +1 one has 

A u 
u 


A W+ 

A l ~ 


- and 


AT 


AT 


A u 
u 

Ad 


( 12 ) 


(13) 


evaluated at x = 0.435. Therefore these measurements will allow a fairly clean flavor 
separation, both for quarks and antiquarks, for some interesting ranges of x— values. We 

2 Actually this theoretical prediction has never been checked experimentally since the single-jet cross 

sections are always given for y = 0 assuming a given pseudo-rapidity cut, usually \y\ < 0.7 or more. 
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see in Fig.7a that A^ + , which is driven by the u and d polarizations, leads to similar 
predictions for all cases. This is mainly due to our knowledge of A u/u, except for x > 0.3 
where it comes out to be larger for BS-h |J. In Fig.7b, for Ad// which is sensitive to the 
d and u polarizations, we see that the various predictions leads to the same general trend, 
with some difference in magnitude due to a large uncertainty in the determination of A d/d. 
Also BS-h |9i] has assumed a larger negative A u/u which is reflected in the behaviour near 
y = —1. The statistical errors have been calculated with a rapidity resolution Ay = 0.2 
and taking into account only the events from the leptonic decay modes. They are smaller 
for W + production which has larger event rates. 


Let us now turn to the parity violating double helicity asymmetries A p )/ and A/// defined 
in eq.(7). Their explicit expressions are given in ref.|23| and we recall that in pp collisions 
they have the following symmetry properties 


A P L I(y)=A p L Y(-y) and A p L v L {y) = -A%(-y). (14) 


A priori these asymmetries are independent observables, but if one makes the reasonable 
assumption AuAd <C ud for all x, one gets the following relations 

A H(y) = A L(y) + A L (-y) , Aj%(y) = A L (y) - A L (-y). (15) 


Therefore these double helicity asymmetries do not contain any additional information, 
but they help to emphasize the differences among the four predictions for A^/ + and A// , 
as seen in Fig.8a,b for A//[ . From an experimental viewpoint, they might be also more 
useful, if one does not have such a small rapidity resolutions. 


Finally let us consider the realistic possibility of having proton-neutron collisions at 
RHIC, which can be achieved with either deuteron or helium-3 beams. The parity violat¬ 
ing single helicity asymmetry Ai for W + production in pn collisions reads 


24 


A L (pn -> W + ]y) 


Au{x a , Mw) u(x b , Mw) - Ad{x a ,M/ v )d{x b ,M/ v ) 
u(x a , Mw) u(x b , M^) + d(x a , M^) d(x b , M^) 


(16) 


assuming the proton is polarized. Of course one has a similar expression for W~ produc¬ 
tion. The double helicity asymmetries are also related to the sum or the difference of two 
Al but this time for A^\ , instead of eq.(15), we have the sum 

AH {Pn —> W + ; y) = A L {pn ^W + ;y) + A L {pn W~ ; -y) (17) 


and therefore the following symmetry property 

AH (Pn W + ; y) = A P L V L {pn - W~ ; -y) (18) 

Similarly, we find for A P Y that eq.(17) holds, with the difference instead of the sum, 
and the symmetry property 

All (m ^W + -y) = - All (pn - W- ; -y) (19) 

Therefore eq.(17) allows us to calculate A p \ , even if one cannot have polarized neu¬ 
trons. We show in Figs.9a,b the four predictions as in the previous cases []. 

3 We realize that since neutron beams are not directly available, but they require deuteron or heliunr-3 
beams, y/s = 500 GeV is obviously too high to be accessible at RHIC. 












6 Concluding remarks 

In this paper we have briefly reviewed the specific properties of some sets of polarized 
parton distributions which have been considered in the recent literature, to describe polar¬ 
ized deep-inelastic scattering data. We have stress the importance of the future polarized 
pp collider at the BNL-RHIC, for studing these distributions, in a kinematic range so far 
inexplored and to provide the first serious tests of the QCD spin sector. 

We have seen that because of the high luminosity of the machine, we will be able to 
pin down the magnitude and the sign of A G(x, Q 2 ) which is, at present, the less well- 
known distribution. This will allow to verify directly one possible interpretation of the 
observed violation of the Ellis-Jaffe sum rules || for proton and neutron, in terms of the 
axial anomaly. For this particular topic, one requires detailed measurements of the double 
hclicity asymmetry All for single jet and direct photon production. These will be fairly 
accurately determined as we have dicussed from the small size of the expected statistical 
errors. 

The copious production of W ± gauge bosons at a/s = 500 GeV will also provide an 
extremely valuable source of information on these distributions and in particular will help 
us to perform a clean flavor separation between u and d quarks and similarly for antiquarks 
u and d. This can be achieved from the measurement of the single helicity asymmetry Al, 
but we have also discussed some interesting features of the double helicity parity violating 
asymmetries A/} and All which are simply related to Al- Although W ± production will 
be easier in pp collisions, we have also considered some interesting aspects of the parity 
violating asymmetries in the plausible situation where we might have pn collisions. 
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Figure captions 


Fig. la The gluon polarization versus x at Q 2 = 4 GeV 2 from ref. || (soft gluon BS-s is 
dotted curve, hard gluon BS-h is dashed curve) and from ref.[|Kj] (BS-i solid curve). 


Fig. l.b The gluon polarization versus x at Q 2 = 4 GeV 2 from ref.JT^] (gluon A GSA is 


solid curve, gluon B GSB is dashed curve, gluon C GSC is dotted curve) and from ref. [|T4 
(standard scenario GRSV, dotted-dashed curve 


Fig. 2 The Q 2 dependence of A^(x,Q 2 ) as predicted from ref. |L4| at fixed values of x. 


Fig. 3.a The double hclicity asymmetry A/ L (p = 0) at p = 0 versus p/ for y/s = 500 
GeV. The solid curves, which are the exact calculations using BS-h (upper curves) and 
GRSV (lower curves), are compared with the dotted curves obtained using eq.(9). 


Fig. 3b Same as Fig. 3a using BS-h (solid curve), BS-s (large dashed curve), BS-i (dotted 
curve), GRSV (dashed-dotted curve) GSA, GSB and GSC (small dashed curves). 

Fig. 4 Various predictions for the p distribution of A V LL at p/ = 40 GeV/c and for 

y/s = 500 GeV (curve labels as in Fig.3b). The statistical errors are calculated for 

A// = 0.5 per bin and for a p/ acceptance A p/ = 10 GeV/c. 

Fig. 5 The double hclicity asymmetry A J //(p = 0) at p = 0 and y/s = 500 GeV versus 
Pjf*, for various parametrizations (curve labels as in Fig.3b). The statistical errors are 
calculated for A?/ = 1 per bin and a p/ 1 acceptance A p/ 1 = 10 GeV/c. 

Fig. 6 Various predictions for the p distribution of A J //(7]) at p = 60 GeV/c and 

y/s = 500 GeV (curve labels as in Fig.3b). The statistical errors are calculated for 

A p = 0.2 and a pi/ 1 acceptance A pi/ 1 = 10 GeV/c. 

Fig. 7a Various predictions for the y distribution of A}// + (y) at y/s = 500 GeV. BS-h 
(solid curve), GRSV (dashed-dotted curve), GSA (dashed curve), GSC (dotted curve). 
The statistical errors are calculated for A p = 0.2. 


Fig. 7b Same as Fig.7a for W production 

Fig. 8a Various predictions for the y distribution of A// (y) for W + production at 
y/s = 500 GeV (curve labels as in Fig.7a). 

Fig. 8b Same as Fig. 8 a for W production. 

Fig. 9a Various predictions for the y distribution of A// (y) for W + production in pn 
collisions at y/s = 500 GeV (curve labels and statistical errors as in Fig.7a). 

Fig. 9b Same as Fig.9a for A/J/y). 
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